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Hitting times for discrete quantum walks on graphs give an average time before the walk reaches 
an ending condition. To be analogous to the hitting time for a classical walk, the quantum hitting 
time must involve repeated measurements as well as unitary evolution. We derive an expression for 
hitting time using superoperators, and numerically evaluate it for the discrete walk on the hypercube. 
The values found are compared to other analogues of hitting time suggested in earlier work. The 
dependence of hitting times on the type of unitary "coin" is examined, and we give an example 
of an initial state and coin which gives an infinite hitting time for a quantum walk. Such infinite 
hitting times require destructive interference, and are not observed classically. Finally, we look at 
distortions of the hypercube, and observe that a loss of symmetry in the hypercube increases the 
hitting time. Symmetry seems to play an important role in both dramatic speed-ups and slow-downs 
of quantum walks. 

PACS numbers: 



I. INTRODUCTION 

Algorithms based on random walks have been im- 
mensely successful in computer science. The most ef- 
ficient solution for 3-SAT is based on the hitting time 
of a classical random walk 

Q- The y 

are also useful in 

algorithms for k-SAT and probability amplification. It 
therefore seems reasonable that their quantum analogues 
might be useful in developing quantum algorithms. This 
has been one of the main motivations for defining and 
analyzing quantum walks. Most quantum algorithms 
to date are based on the Quantum Fourier Transform 
(QFT), like Shor's algorithms for factoring and discrete 
log 2J. But the use of QFT seems to be restricted. For 
example, it does not seem to be useful in some non- 
Abelian hidden subgroup problems like the Graph Iso- 
morphism problem. Therefore, there is a need for new al- 
gorithmic tools for the analysis of such problems. This is 
another reason to study the properties of quantum walks 
on graphs. 

Quantum walks come in two versions: continuous time 
and discrete time. Continuous time versions have been 
studied in 

0- H, II- In H, it has bee n demonstrated 
that for the "glued trees" graph, the quantum walk is ex- 
ponentially faster than a classical random walk. Discrete 
time versions have been studied extensively for the walk 
on the line @, 0, 0, and the hypercube El, [3, 
|l3|. Quantum walks for any general (irregular) graphs 
have been defined in [l^ and In it was shown 
that there is a quadratic speed up with discrete quan- 
tum walks in searching an unstructured database for a 
marked node. [17] describes how these concepts can be 
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used for the problem of element distinctness. Reviews on 
the recent developments in this field can be found in [lfij 
and 0. 

Classical random walks are characterized by a number 
of time scales: mixing time, cover time, correlation time, 
and hitting time. Similar quantities can be defined for 
quantum walks, analogous to their classical counterparts, 
though some types of classical behavior are not guaran- 
teed to hold for quantum walks. In El a lower bound on 
the mixing time for quantum walks was found in terms 
of the conductance of the underlying gra ph. Some no- 
tions of hitting time have been defined in [Tj , where an 
upper bound on the hitting time has been found for the 
quantum walk on the hypercube. 

In a classical discrete random walk, a hitting time is 
the average time for a walk beginning on a particular 
starting vertex to arrive for the first time at a particular 
ending vertex or group of vertices. It is not sufficient 
for the ending vertex merely to have positive probability: 
the particle must actually arrive there. Therefore, hitting 
times are generally longer than the shortest path between 
the starting and ending vertices. An analogous definition 
can be made for a continuous random walk, but in this 
paper we will only consider discrete walks. 

This definition becomes complicated for a quantum 
walk, where the particle can be in a superposition of 
many vertex locations at once. To capture the idea of 
first arrival time, our definition of hitting time is the ex- 
pected time for the particle to be at the final vertex for 
the first time as determined by a measurement of the 
final vertex. At every time step we perform a measure- 
ment to see if the particle is at the final vertex or not. 
This notion of hitting time seems to be the most natural 
extension of the definition of hitting time to the quantum 
case. Using this definition, we will obtain an expression 
for hitting time on any graph, and produce numerical re- 
sults showing that it is orders of magnitude lower than 
the classical hitting time for the quantum walk on the 
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hypercube. 

This paper is organized as follows. In section II, a 
quantum walk on a regular graph is defined and this def- 
inition applied to a walk on the hypercube. In section III, 
we review the classical definition of hitting time. Then 
we define the quantum hitting time and derive an ex- 
pression for it on any regular graph (or on an irregular 
graph when the walk is appropriately defined). We also 
briefly review other definitions of hitting time that have 
been used for quantum walks. In section IV, we discuss 
the results of simulations for a walk on the hypercube 
using the Grover coin, and compare them to the classi- 
cal hitting time and the other definitions for a quantum 
hitting time, showing that in this case the hitting time is 
far shorter than these other definitions. Then in section 
V we show that the hitting time for a quantum walk can 
be much longer than the classical hitting time for certain 
combinations of initial state and coin, and can even be 
infinite. We demonstrate this for a walk on the hyper- 
cube with the discrete Fourier transform (DFT) coin. In 
section VI we look at a distorted version of the hyper- 
cube, and show that the loss of symmetry increases the 
hitting time. We conclude in section VII with a discus- 
sion of the role of symmetry in the speed-up of quantum 
walks. 



II. QUANTUM WALK ON THE HYPERCUBE 
A. Quantum walks 

A discrete time quantum walk is defined as the appli- 
cation of a unitary operator [Zona Hilbert space repre- 
senting the graph; each vertex of the graph has a set of 
basis states associated with it, and the unitary operator 
can only cause a transition between two basis states if the 
vertices associated with them are connected by an edge. 
If the graph is regular, the Hilbert space can be a tensor 
product of two Hilbert spaces: in a walk on a regular 
graph with N vertices and degree d, U is applied to an 
iVd-dimensional Hilbert space W W c , where W is the 
Hilbert space of the position (i.e., the vertex), and H c is 
the Hilbert space of a G?-sided "coin," which is "flipped" 
at each step to determine which edge to walk on. 

For a given graph, if we label each edge incident on a 
vertex by a number from 1 through d, such that an edge 
connecting two vertices gets the same label at both ends, 
then each vertex is a basis state in the position space and 
each direction is a basis state in the coin space. There- 
fore, the basis states for the position and coin spaces can 
be labelled |0), |2), . . . ,\N - 1) and |0), |2), . . . , \d - 1). 
The unitary U is of the form 



U = SF, 



(1) 



where we call S the shift matrix and F the coin flip ma- 
trix. S is a permutation matrix which shifts the parti- 
cle from its present vertex along the edge indicated by 



the coin state, in a way analogous to a classical ran- 
dom walk. The coin flip matrix F acts solely on the 
coin space, so it is of the form I (g> C. C can, in prin- 
ciple, be any unitary matrix, though in general we con- 
centrate on examples with some kind of useful structure. 
The most common coins analyzed are the Grover coin G, 
the Discrete Fourier Transform (DFT) coin D, and the 
Hadamard coin H. While the Grover and the DFT coins 
exist for all dimensions, the Hadamard coin exists only 
for dimensions 2™ for some n. These coins are given by 
the following matrices: 



6 = 2|tf)<*| -I, 
where |*> = ^ES), i.e., 

G = 



(2) 



2 



1 \ 

i 

d 



D = 



1 

7d 



V 



A i 

1 LJ 



(3) 



1-V 



1 \ 

,d-i 



\1 c^- 1 u; 2 ^- 1 ) 
where lo — exp(27ri/c?), and 

where 



(4) 



H: 



2. 



V2 V 1 - 1 



(5) 



(6) 



The shift operator is applied after the coin operator. The 
shift operator moves the particle from a vertex along the 
edge given by the direction number of the coin state i.e., 



(7) 



where v(i) is the vertex connected to v via the edge num- 
bered i. It is interesting to note that the term "quantum 
random walk" is, in a sense, a misnomer, because the 
randomness in a quantum walk is introduced by quantum 
measurements where one of the measurement outcomes 
takes place at random. Thus, there is no extra random- 
ness in the walk, i.e., outside of what is introduced by 
quantum mechanics itself. 



B. The Hypercube 

The hypercube of dimension n is a set of 2™ vertices 
each with a degree n. These vertices can be numbered 
by a n-bit string from (00 ... 0) through (11 . . . 1). Each 
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edge leads in a particular direction, and can be labeled 
by an integer j, < j < n. Adjacent vertices are the ones 
whose bit assignments differ by a single bit. For example, 
for the hypercube in 3 dimensions, the vertices (Oil) and 
(111) are adjacent to each other, connected by an edge in 
the 2 direction. We can treat the bit strings of the ver- 
tices as n-dimensional Boolean vectors, whose elements 
are all or 1 ; the directions of the edges then correspond 
to vectors ej with element j equal to 1 and all other ele- 
ments zero. The vertex adjacent to a given vertex \x) in 
the j direction can be labeled \x © e~l) ie{0, 1, . . . , n — 1}. 

Since the vertices of the hypercube can be labeled |0) 
through 1 2" — 1) and the directions labeled |0) through 
|n-l), the Hilbert space is H 2 " ® The shift operator 
is of the form 



EE 

j=0 x 



x®e-j,j)(x,j\, 



(8) 



where \v, j) — \v) (g) \ j), \v) is the position state and \j) 
the coin state. 



time starting at any vertex with hamming weight x, then 

77 — T T. 

t(x) = t(x + 1) + -t(x - 1) + 1, (10) 

n n 

with the boundary condition r(n) = 0. This simplifies to 



„/n n — x — 1 . , \ n 
A(x) = _A(x + l) 



x + 1 



x + l J 



(11) 



where A(x) = t(x) — t(x + 1). Using this recursive for- 
mula, we obtain 



n-l 



r ( 0) = ^A ( x) = ^ ^„!7 + , (12, 



x=0 



x=0 



C: 



This sum can readily be evaluated for reasonable sizes of 
n; we use this expression to compare the classical hitting 
time to the quantum hitting time. We define the hitting 
time of a quantum walk next. 



B. Quantum hitting time 



III. HITTING TIME 

A. Classical hitting time 

Given a regular undirected graph and a particle which 
starts at some vertex, the classical random walk is de- 
fined as follows. At each vertex, the particle moves along 
any edge incident on the vertex with some predefined 
probability. This procedure is then repeated at the new 
vertex. The walk continues until the particle arrives at 
("hits") a certain vertex (called the "final vertex") for 
the first time. The hitting time is defined as the average 
time until the particle hits the final vertex: 

oo 

T(v)=J2tPv(t), (9) 
t=0 

where t(v) is the hitting time given that the walk starts 
at vertex v and p v (t) is the probability that the particle 
hits the final vertex for the first time at time step t (first 
crossing probability) given that it was at v at t = 0. 

Let us now specialize to the case of the hypercube, 
where the the final vertex is assumed to be 11 ■ ■ ■ 1. We 
would like to find the hitting time starting from 00 • • • 0. 
For the classical walk on the hypercube, one can arrive at 
a recursive relation involving the hitting time. First, from 
the symmetry of the hypercube one can conclude that 
the hitting time depends only on the hamming weight 
of the starting vertex rather than the vertex itself. The 
hamming weight is the number of l's in the string of 
bits. At hamming weight x, there are C" = n\/x\(n — x)\ 
vertices. The probability to walk to a vertex with weight 
x + 1 is (n — x)/n, and the probability to walk to a vertex 
with weight x — 1 is x/n. So, if t(x) denotes the hitting 



We define the hitting time of a quantum walk in close 
analogy to that of a classical random walk. For a classical 
walk, hitting time is the average time taken for the parti- 
cle to hit the final vertex for the first time. To carry this 
over to the quantum case, we must give a proper meaning 
to the phrase "for the first time." The only reasonable 
way to do this is to measure the final vertex at every step 
to see if the particle has arrived or not. Without such 
a measurement at every step, one cannot appropriately 
define the first crossing probability. To define the hitting 
time this way, we first define what is called a measured 
quantum walk on a graph. 

Suppose we have a walk on a regular graph of degree d 
with N vertices, where the particle begins at a vertex xo 
and walks till it reaches the vertex x / . Assume that the 
coin starts in the initial state |co). The initial state of the 
particle is then,|*) = \x ) <g> |c ); po = I^X^I where po 
is the density operator corresponding to the initial state 
\^). The measured walk is now the alternating applica- 
tion of a unitary evolution operator U (the product of the 
shift and the coin flip operators) and a projective mea- 
surement M with 2 outcomes Pf and Qf = I — Pf, where 
Pf = \xf)(xf\ ® I is the projector onto the final vertex 
for any coin state. Thus, if the particle is not found in 
the final vertex after t time steps, the state is: 

= (QfUMiQfun (13) 

11 Tr{(Q f Uy P ((Q f UnY 
If the particle is found in the final vertex, the walk is 
assumed to end. The state \xf) becomes an absorbing 
boundary for this (measured) walk. Now, the first cross- 
ing probability at time step t can be defined (analogous 
to the classical case) to be the following expression: 

p(t) = TriPfUlQf^polWQfY^mPf}. (14) 
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With this definition for the probability, the hitting time 
becomes 



(15) 



To sum this series requires a slightly different expression 
for the first arrival probability pit). We rewrite (|14|1 in 
terms of the following super operators (i.e., linear trans- 
formations on operators): 



Np = QfUpWQ f 

y P = PfUpWP f . 



(16) 



In terms of TV and y, pit) = Tr{yN t ^ 1 po}- We now 
evaluate the hitting time by introducing a new superop- 
erator, 



(17) 



which is a function of a parameter I. The hitting time 
now becomes 



dl 



Tr{yOil)p } 



(18) 



If the superoperator X — IN is invertible, then we can 
replace the sum (I17|) with the closed form 



o(i) = iii - iAfy 



(19) 



If X — lj\f is not invertible, then this sum instead is a 
pseudoinverse — that is, the inverse restricted to the sup- 
port of the superoperator. So long as the graph in ques- 
tion is finite, this is well-defined. The derivative in ltTH|) 
is then 

^(1) = il-U^+Mil-M)- 2 = iX-Af)- 2 . (20) 
This gives us the following expression for the hitting time: 



T = TY{yiX-Af)- 2 p }. 



(21) 



The meaning for the terms in the above expression (e.g. 
X — J\f)~ 2 ) can be given by vectorizing all the operators. 
Any matrix can be vectorized by turning its rows into 
columns and stacking them up one by one, so that a 
D x D matrix becomes a column vector of size D 2 . For 
example: 

/an\ 
ai3 

'an ai2 ai 3 \ a 2 i 

^21 «22 123 — ► «22 
V a 31 a 32 033/ «23 
031 
032 
V a 33/ 



Consequently the superoperators become matrices of size 
D 2 x D 2 . 

Denoting the vectorized quantities as p v , etc., we ob- 
tain the hitting time as 



T = Tr v {y(I-M)- 2 p v }, 



(22) 



where {ypf = iP f UpWP f ) v and iAfp) v = 

iQ f UpU^Qf) v and Tr u {.} is the equivalent of the trace 
operation for vectorized quantities. It is just the inner 
product of the resultant vector with the vectorized 
identity matrix: Tr v {O v } = T° ■ O v . 

[ljj, we obtain 



Using Roth's lemma 

iNpY = [(QjU)®iQ f Uy 

(yp) v = \(PfU)®{PfU)' 



p 

p v - 



(23) 



Therefore, the expression for hitting time becomes 



F ■ (iP f tJ) % iP f UriX-iQ f U) ® (<3/?7)*)~V 



(24) 



C. Other definitions of hitting time 

Two definitions of the quantum hitting time were given 
111 one-shot hitting time and concurrent hitting time. 
The one-shot hitting time is defined for an unmeasured 
walk. It is the time at which the probability of being in 
the final state is greater than some given value. More 
precisely, given some probability p, the one-shot hitting 
time is defined as the lowest time T such that 



\(x f \U T \x )\ 2 >p 



(25) 



where xj and xq are the final and initial states, and U 
is the evolution operator (as defined above). Essentially 
this same definition of hitting time was used in the analy- 
sis of the continuous-time walk on the hypercube in [2(J . 
This definition is useful if it is known that at some time 
the probability to be in the final state will be higher than 
some reasonable value; but for a general graph, this is not 
guaranteed. 

The concurrent hitting time, by contrast, is defined 
for a measured quantum walk. Given a probability p, 
it is the the time T such that the measured walk has a 
probability greater than p of stopping at a time less than 
T . It has been proved that the concurrent hitting time 
is T = -| for p — &X n ; gii n ) for a hypercube of dimension 
n. Since we consider only the measured quantum walk 
in this paper, we compare our numerical results to the 
numerical simulation of the concurrent hitting time and 
the bound on it derived in In the next section, this 
is redefined in terms of the residual probability 1 — p and 
plotted against the hitting time defined in the previous 
section. 
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Dimension of hypercube 



FIG. 1: Hitting times of Classical and Quantum walks 
(semi-log scale) 



FIG. 3: Comparison of Quantum hitting time, stop- 
ping time and the bound on the stopping time for 
dimensions 50 to 60(log scale) 



If we think of quantum walks as a possible route to new 
algorithms, then the concurrent hitting time corresponds 
to the time needed to find a solution with probability 
greater than p. The definition of hitting time used in 
this paper corresponds more to a typical running time 
for the algorithm. Both definitions could prove useful for 
particular purposes. 



IV. NUMERICAL RESULTS 
A. Results for the Grover coin 



10 s 
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We calculated the hitting time by evaluating the above 
expression H21(l in Matlab, and comparing the result to 
a lower bound obtained by iterating the quantum walk 
for a large number of steps. Because of the multiple 
tensor products, the size of the matrices Af and y is 
(2™n) 2 x (2 n n) 2 . 

We can make this computation more tractable in the 
case where the coin-flip unitary is the Grover coin G 
defined in (J3J, and the specific starting state is p — 
|0)(0|(g>|^)(</>|, where \<j>) = -^X)2=il*)- It was shown by 
Shenvi et al. 0] that for the walk on the hypercube with 
this initial state and the Grover coin, the state remains 
always in a 2n-dimensional subspace, where the walk is 
on a line with n + 1 points. This is rather similar to the 
simplification made in the classical case, when we kept 
track only of the Hamming weight of the current vertex. 
With this simplification, the operators in l121[| reduce to 
(2n) 2 x (2n) 2 matrices, which makes explicit calculations 
possible even for high dimensional hypercubes. 

We write a set of basis states for this subspace as 
|-R,0},|L,l),|.R,l},...,|.R,n - l),\L,n), where the first 
label says whether the state is "right-going" or "left- 
going," and the second label gives the Hamming weight 
of the state. The initial state is \R, 0), and the final state 
is \L,n). (Note that there are no states |£,0) or \R,n).) 
Restricted to this subspace, the S and C matrices become 

TO 

S = ^2 \R,x){L,x + l\ + \L,x+ 1) (R, x\ (26) 

x=0 



FIG. 2: Comparison of Quantum hitting time, stop- 
ping time and the bound on the stopping time for 
dimensions 10 to 20(log scale) 



C\L,x) = — coslo x \L, x) + sinu) x \R, x), 
C\R, x) = sinw x \L, x) + cosuj x \R, x), (27) 

where cosu; x = 1 — 2x/n. We see that for the walk in 
this subspace, the coin flip is no longer independent of 
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the position; this is quite analogous to the reduction of 
the classical walk to the Hamming weight, in which the 
probabilities favor walking towards x = n/2. 

Now, going back to the definition of the hitting time 
as a sum of series in Eq. (|15|) where the probability is 
given by Eq. (|14J) . we can define the e-stopping time 
(r e ) as the time step T at which the total probability 
Y^J—i > 1 — £ i where p(i) is as in Eq. I|14|). We calcu- 
late an estimate of the hitting time (r est (e) as a function 
of e by summing the series in Eq. i|15[) up to r e . Figure ^ 
shows the classical and quantum walks on the hypercube 
for 100 dimensions. The exact hitting time r calculated 
by computing the expression in equation 1|24|) is plotted as 
the dotted line and T es t(e) for e = 0.001 is plotted as the 
solid line. These two lines almost coincide in the graph. 
The hitting time on the Y-axis is plotted in log scale. We 
can see that hitting time for the quantum walk is a low 
order polynomial whereas the classical walk is exponen- 
tial. There is a very dramatic speed-up in the quantum 
case. 

Figure|21plots the quantum hitting time, the e-stopping 
time and the bound on the e-stopping time (obtained 
in against e for dimensions from 10 to 20. Both 

the axes are in log scale. Figure |3 plots the same for 
dimensions from 50 to 60. It can be seen from these 
two figures that the r e and the bounds on it become less 
tight for higher dimensions. Clearly, the average walk 
ends much faster than these bounds might suggest. 

B. Results for the DFT coin 

The hitting time for the hypercube using the DFT coin, 
by contrast to the Grover coin case considered above, can 
actually be infinite. For n=4 we will demonstrate that 
for the same initial condition as with the Grover coin, 
the hitting time for a quantum walk using the DFT coin 
is infinity. This is because there exist eigenvalues of the 
evolution operator hi whose eigenvectors have an overlap 
with the initial state, but have no overlap with the final 
vertex for any state of the coin. 

Suppose there are states \<j)) which have no overlap with 
the final vertex, (<f)\Pf\(f)) = 0, and which are eigenstates 
of the evolution operator U : U\<j>) — exp(i8)\<j>) . If the 
system is in the state \<f>), clearly there is no probability 
to ever detect the particle in the final vertex. Let V 
be a projector onto all such states \<f>). Then VPf — 
PfV — 0, and [P,U] = 0. One can write the initial 
state as a superposition of vectors in this subspace and 
its orthogonal complement: 

|*) = P|#) + (±-P)|#). (28) 

Any state that begins in the subspace selected by V will 
remain there for all time, and any state in the orthog- 
onal complement will stay there; this follows from the 
fact that the projectors commute with both the unitary 
transformation U and the measurement operator Pf. As 



one starts the walk, the probability that the particle 

- t „ » 

never reaches the final state is (^|f/T pjj * | \|7) , which is 
(#|P|¥>. 

In order for this probability to be nonzero, there must 
be eigenstates of the unitary evolution operator U which 
have no amplitude for the final vertex. We can readily 
demonstrate this for the hypercube with the DFT coin. 
Consider the 4-dimensional hypercube. Numerically di- 
agonalizing the evolution operator U given by we find 
it has i, —i, 1 and —1 among its eigenvalues, each with a 
degeneracy of 8. Since the subspace corresponding to 
the final vertex is 4-dimensional, it is clearly possible to 
construct a superposition of eigenvectors of any of these 
eigenvalues so that it has no overlap with the final ver- 
tex in any coin state. For each of the four degenerate 
eigenvalues we can construct a 4-dimensional subspace 
of eigenvectors with no overlap with the final vertex, giv- 
ing a 16-dimensional space for all such eigenvectors. By 
numerically constructing an orthonormal basis for this 
space, we can find an expression for the projector V and 
measure its overlap with the initial state. 

We considered in particular the initial state where the 
particle was located at the |00 . . . 0) vertex and the coin 
is in an equal superposition of basis states .. . , \n). 
For the hypercube with n — 4 and the given initial state, 
the probability (\I r |'P|\Ij) is 0.4286, which exactly matches 
the total probability to never hit the final node after a 
large number of iterations in our numerical simulations. 
Thus, the probability is close to half that the particle 
never reaches the final state and the hitting time becomes 
infinity. 

This demonstrates a property of quantum walks not 
seen in their classical counterparts: for certain initial 
conditions, there is a nonzero probability that the parti- 
cle never reaches the final state, even though the initial 
and final states of the graph are connected. For a quan- 
tum walk with substantial degeneracy, this phenomenon 
is likely to be generic. It might be possible to make 
the hitting time finite by choosing an appropriate initial 
condition — clearly this happens for the Grover coin — but 
for some coins this may require an initial condition which 
is not localized on one vertex. From our simulations, it it 
seems that for higher dimensions the DFT coin behaves 
similarly to n = 4. For example, for n = 5, our simula- 
tions show that the probability to hit the final node in- 
creases slowly but does not reach 1 even after many time 
steps. This could be due to the fact that the final vertex 
has no overlap with some eigenvectors of the evolution 
operator (as for n = 4), and additionally that it overlaps 
very little for some other eigenvector. This would make 
the probability increase slowly but never reach 1. 

C. Results for a quantum walk on a distorted 
hypercube 

If, as seems likely, the dramatic speed-ups (and slow 
downs) of quantum walks over their classical counter- 
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FIG. 4: Comparison of hitting times on the regular 
and distorted hypercubes 

parts depend on the symmetry of the graph, it should be 
instructive to see the effect of deviations from that sym- 
metry. In this section, we look at results for the mea- 
sured walk using the Grover coin on a distorted hyper- 
cube. The distorted hypercube is defined by construct- 
ing the usual hypercube, and then switching two of the 
connections. Pick 4 vertices which form a face-for exam- 
ple, (0 . . . 00), (0 . . . 01), (0 . . . 10), (0 ... 11). Calling these 
vertices A, B, C, D for short, we distort the hypercube by 
connecting A to D and B to C, and removing the edges 
between A and B and between C and D. This is still a 
regular graph, and the same quantum walk can be used 
without having to redefine the evolution operator. Unlike 
the usual hypercube, it is no longer a bipartite graph, and 
the walk can no longer be reduced to a walk in Hamming 
weight. 

Figurc0]plots the hitting time of a quantum walk on a 
distorted hypercube together with that of a classical walk 
and a quantum walk on regular hypercubes for compari- 
son. The hitting time for a quantum walk on a distorted 
hypercube is more than that of a quantum walk on a 
regular hypercube, but still is much smaller than the hit- 
ting time of a classical walk. In fact, as the dimension 
increases one can see that the hitting times of the quan- 
tum walk on the distorted and regular hypercubes con- 
verge towards each other. This presumably reflects the 
fact that for higher dimensions the symmetry is mostly 
unchanged. 

V. CONCLUSIONS 

We have examined the definition of hitting time for 
measured quantum walks and analyzed some of its prop- 
erties. We used this definition to obtain an expression 
for hitting time which is valid on any general graph as 
long as the unitary evolution operator U of the walk is 



defined. We simulated this hitting time for a measured 
quantum walk using the Grover coin and compared it 
to the classical hitting time and to the bounds obtained 
on it; the quantum hitting time is exponentially smaller 
than the classical hitting time. We also showed that the 
bounds on the hitting time obtained in ^lj become less 
tight as the dimension increases. 

However, simply making a walk quantum does not 
guarantee a speed-up over the classical case. We demon- 
strated that the hitting time for quantum walks can de- 
pend sensitively on the initial condition, unlike classical 
walks. For certain initial states, the DFT walk can have 
infinite hitting time, a phenomenon not possible in clas- 
sical random walks. This dependence on the initial state 
varies with the coin used, since for the same initial state 
the Grover walk has a polynomial hitting time. This in- 
finite hitting time is directly related to the degeneracy of 
the eigenvalues of the evolution operator. If the evolu- 
tion operator is highly degenerate, then it is very likely 
that there exist initial states which give infinite hitting 
times. 

While the exact cause of the speed-up in quantum hit- 
ting time is not completely clear, it seems very likely 
that the symmetry of the graph plays a major role in 
both the speed-up and slow down of the quantum walk. 
In the faster quantum walk, the different paths leading 
to the final vertex interfere constructively, enhancing the 
probability of arrival; paths which lead to "wrong" ver- 
tices interfere destructively, reducing the probability of 
meandering around in the graph for long times. Unlike a 
classical random walk, the quantum walk is sensitive to 
the presence of a global symmetry which is not apparent 
at a purely local level. This phenomenon probably leads 
to the speed-up of the continuous-time quantum walk on 
the glued-trees graph as well 

However, this same symmetry is undoubtedly the cul- 
prit in the slow-down observed for the DFT walk. The 
existence of states which never arrive at the final vertex 
is made possible by the degeneracy of the evolution oper- 
ator U — a degeneracy which arises due to the symmetry 
of the graph. The existence of states which never arrive 
at the final vertex can also be seen as an interference ef- 
fect, only in this case the interference of paths which lead 
to the final vertex is destructive: all amplitude to make 
a transition to the final vertex cancels out. 

This interpretation is supported by the quantum walk 
on the distorted hypercube. We observe that the hitting 
time is worse than that of the usual hypercube, but still 
much smaller than that of a classical walk. The curve 
of the hitting time on the distorted hypercube seems to 
converge slowly to that of a quantum walk on the reg- 
ular hypercube. This is probably because the distortion 
(described in section 4) is very mild. As the dimension 
grows, and with it the number of edges and vertices, this 
distortion has less effect on the overall symmetry. 

We hasten to add that symmetry of the graph is not 
the sole reason for speed-ups in quantum walks. A poly- 
nomial speed-up has been demonstrated in the quantum 
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walk versions of the search of an unstructured database 
[l2| and the element distinctness problem [17|. However, 
the dramatic exponential speed-ups have all been demon- 
strated in highly symmetric graphs. 

This suggests that the most promising direction to look 
for new algorithms based on quantum walks is for prob- 



lems which possess a global symmetry, but for which this 
symmetry is not apparent at the level of individual can- 
didate solutions. Yet even for such a problem, care will 
have to be taken if quantum mechanics is to serve as a 
blessing and not as a curse. 
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